Ancestry dynamics and trait selection in a designer cat breed.

The Bengal cat breed was developed from intercrosses between the Asian leopard cat, Prionailurus bengalensis, and the domestic cat, Felis catus, with a last common ancestor approximately 6 million years ago. Predicted to derive ∼94% of their genome from domestic cats, regions of the leopard cat genome are thought to account for the unique pelage traits and ornate color patterns of the Bengal breed, which are similar to those of ocelots and jaguars. We explore ancestry distribution and selection signatures in the Bengal breed by using reduced representation and whole-genome sequencing from 947 cats. The mean proportion of leopard cat DNA in the Bengal breed is 3.48%, lower than predicted from breed history, and is broadly distributed, covering 93% of the Bengal genome. Overall, leopard cat introgressions do not show strong signatures of selection across the Bengal breed. However, two popular color traits in Bengal cats, charcoal and pheomelanin intensity, are explained by selection of leopard cat genes whose expression is reduced in a domestic cat background, consistent with genetic incompatibility resulting from hybridization. We characterize several selective sweeps in the Bengal genome that harbor candidate genes for pelage and color pattern and that are associated with domestic, rather than leopard, cat haplotypes. We identify the molecular and phenotypic basis of one selective sweep as reduced expression of the Fgfr2 gene, which underlies glitter, a trait desired by breeders that affects hair texture and light reflectivity.

Independent COL17A1 Variants in Cats with Junctional Epidermolysis Bullosa.

Epidermolysis bullosa (EB), characterized by defective adhesion of the epidermis to the dermis, is a heterogeneous disease with many subtypes in human patients and domestic animals. We investigated two unrelated cats with recurring erosions and ulcers on ear pinnae, oral mucosa, and paw pads that were suggestive of EB. Histopathology confirmed the diagnosis of EB in both cats. Case 1 was severe and had to be euthanized at 5 months of age. Case 2 had a milder course and was alive at 11 years of age at the time of writing. Whole genome sequencing of both affected cats revealed independent homozygous variants in COL17A1 encoding the collagen type XVII alpha 1 chain. Loss of function variants in COL17A1 lead to junctional epidermolysis bullosa (JEB) in human patients. The identified splice site variant in case 1, c.3019+1del, was predicted to lead to a complete deficiency in collagen type XVII. Case 2 had a splice region variant, c.769+5G>A. Assessment of the functional impact of this variant on the transcript level demonstrated partial aberrant splicing with residual expression of wildtype transcript. Thus, the molecular analyses provided a plausible explanation of the difference in clinical severity between the two cases and allowed the refinement of the diagnosis in the affected cats to JEB. This study highlights the complexity of EB in animals and contributes to a better understanding of the genotype-phenotype correlation in COL17A1-related JEB.

Identification of an ADAMTS2 frameshift variant in a cat family with Ehlers-Danlos syndrome.

We investigated 4 European domestic shorthair kittens with skin lesions consistent with the dermatosparaxis type of the Ehlers-Danlos syndrome, a connective tissue disorder. The kittens were sired by the same tomcat but were born by 3 different mothers. The kittens had easily torn skin resulting in nonhealing skin wounds. Both clinically and histologically, the skin showed thin epidermis in addition to inflammatory changes. Changes in collagen fibers were visible in electron micrographs. The complete genome of an affected kitten was sequenced. A one base pair duplication leading to a frameshift in the candidate gene ADAMTS2 was identified, p.(Ser235fs*3). All 4 affected cats carried the frameshift duplication in a homozygous state. Genotypes at this variant showed perfect cosegregation with the autosomal recessive Ehlers-Danlos syndrome phenotype in the available family. The mutant allele did not occur in 48 unrelated control cats. ADAMTS2 loss-of-function variants cause autosomal recessive forms of Ehlers-Danlos syndrome in humans, mice, dogs, cattle, and sheep. The available evidence from our investigation together with the functional knowledge on ADAMTS2 in other species allows to classify the identified ADAMTS2 variant as pathogenic and most likely causative variant for the observed Ehlers-Danlos syndrome.

SOAT1 missense variant in two cats with sebaceous gland dysplasia.

Spontaneously arisen hereditary diseases in domestic animals provide an excellent opportunity to study the physiological functions of the altered genes. We investigated two 4-month-old sibling domestic short haired kittens with dry dark debris around the eyes, nose, and ears, dark crusting on the legs and a thin poor hair coat. Skin biopsies revealed abnormal sebaceous gland morphology with lack of normal sebocyte arrangement and differentiation. Hair follicles had a distorted silhouette, interpreted as a change secondary to the observed sebaceous gland dysplasia. Whole genome sequencing on both affected kittens and 65 genetically diverse feline genomes was performed. Filtering for variants that were present in both kittens but absent from the control genomes revealed a homozygous missense variant in SOAT1, encoding sterol O-acyltransferase 1. The protein is localized in the endoplasmic reticulum and catalyzes the formation of cholesteryl esters, an essential component of sebum and meibum. The identified SOAT1:c.1531G > A variant is predicted to change a highly conserved glycine residue within the last transmembrane domain of SOAT1, p.Gly511Arg. In mice, variants in Soat1 or complete knockout of the gene lead to the "hair interior defect" (hid) or abnormal Meibomian glands, respectively. SOAT1:c.1531G > A represents a plausible candidate variant for the observed sebaceous gland dysplasia in both kittens of this study. The variant was not present in 10 additional cats with a similar clinical and histopathological phenotype suggesting genetic heterogeneity. SOAT1 variants should be considered as potential cause in hereditary sebaceous gland dysplasias of humans and domestic animals.

SDR9C7 missense variant in a Chihuahua with non-epidermolytic ichthyosis.

Ichthyoses represent a heterogeneous group of cornification disorders that are associated with skin barrier defects. We investigated a 9-month-old Chihuahua showing excessive scale formation. Clinical and histopathological examinations revealed non-epidermolytic ichthyosis and a genetic defect was suspected. We therefore sequenced the genome of the affected dog and compared the data with 564 genetically diverse control genomes. Filtering for private variants identified a homozygous missense variant in SDR9C7, c.454C>T or p.(Arg152Trp). SDR9C7 is a known candidate gene for ichthyosis in humans and encodes the short-chain dehydrogenase/reductase family 9C member 7. The enzyme is involved in the production of a functional corneocyte lipid envelope (CLE), a crucial component of the epidermal barrier. Pathogenic variants in SDR9C7 have been described in human patients with autosomal recessive ichthyosis. We assume that the identified missense variant in the affected Chihuahua of this study impairs the normal enzymatic activity of SDR9C7 and thus prevents the formation of a functioning CLE, resulting in a defective skin barrier. To the best of our knowledge, this is the first report of a spontaneous SDR9C7 variant in domestic animals.

Heterozygous ATP2A2 missense variant identified in a Shih Tzu with Darier disease.

Darier disease is caused by heterozygous loss of function variants in the ATP2A2 gene encoding the endoplasmic/sarcoplasmic reticulum Ca2+ pump ATP2A2. Defective intracellular calcium signaling in the epidermis results in a loss of desmosomal adhesion and the development of characteristic skin lesions. In this study, we investigated a Shih Tzu that developed erythematous papules on the ventrum and, over time, the dorsal neck and a nodule in the right ear canal with secondary ear infection. Histopathologic examination demonstrated discrete foci of acantholysis affecting suprabasal layers of the epidermis. Whole genome sequencing of the affected dog identified a heterozygous missense variant, p.N809H, affecting an evolutionarily conserved amino acid residue of the ATP2A2 protein. The highly characteristic clinical and histopathologic findings together with a plausible variant in the only known functional candidate gene establish the diagnosis of canine Darier disease in the studied dog and highlight the potential of genetic analyses as complementary diagnostic approach in veterinary medicine.

A de novo variant in the keratin 1 gene (KRT1) in a Chinese shar-pei dog with severe congenital cornification disorder and non-epidermolytic ichthyosis.

A 3-months old Chinese shar-pei puppy with ichthyosis was investigated. The dog showed generalized scaling, alopecia and footpad lesions. Histopathological examinations demonstrated a non-epidermolytic hyperkeratosis. The parents of the affected puppy did not show any skin lesions. A trio whole genome sequencing analysis identified a heterozygous de novo 3 bp deletion in the KRT1 gene in the affected dog. This variant, NM_001003392.1:c.567_569del, is predicted to delete a single asparagine from the conserved coil 1A motif within the rod domain of KRT1, NP_001003392.1:p.(Asn190del). Immunohistochemistry demonstrated normal levels of KRT1 expression in the epidermis and follicular epithelia. This might indicate that the variant possibly interferes with keratin dimerization or another function of KRT1. Missense variants affecting the homologous asparagine residue of the human KRT1 cause epidermolytic hyperkeratosis. Histologically, the investigated Chinese shar-pei showed a non-epidermolytic ichthyosis. The finding of a de novo variant in an excellent functional candidate gene strongly suggests that KRT1:p.Asn190del caused the ichthyosis phenotype in the affected Chinese shar-pei. To the best of our knowledge, this is the first description of a KRT1-related non-epidermolytic ichthyosis in domestic animals.

KRT5 missense variant in a Cardigan Welsh Corgi with epidermolysis bullosa simplex.

Epidermolysis bullosa (EB) is a group of blistering disorders that includes several subtypes, classified according to their level of cleavage. Typical clinical signs are blisters and erosions resulting from minimal trauma. The disease has been described in many mammalian species and pathogenic variants in at least 18 different genes have been identified. In the present study, we investigated a Cardigan Welsh Corgi with congenital clinical signs consistent with epidermolysis bullosa. The puppy had blisters and erosions on the paw pads, and the oral mucosa. Histologic examination demonstrated the typical clefting between the dermis and epidermis and confirmed the clinical suspicion. We obtained whole genome sequencing data from the affected puppy and searched for variants in candidate genes known to cause EB. This revealed a heterozygous missense variant, KRT5:p.(E476K), affecting the highly conserved KLLEGE motif of keratin 5. The mutant allele in the affected puppy arose owing to a de novo mutation event as it was absent from both unaffected parents. Knowledge of the functional impact of KRT5 variants in other species together with the demonstration of the de novo mutation event establishes KRT5:p.(E476K) as causative variant for the observed EBS.

Independent COL5A1 Variants in Cats with Ehlers-Danlos Syndrome.

We investigated four cats with similar clinical skin-related signs strongly suggestive of Ehlers-Danlos syndrome (EDS). Cases no. 1 and 4 were unrelated and the remaining two cases, no. 2 and 3, were reportedly siblings. Histopathological changes were characterized by severely altered dermal collagen fibers. Transmission electron microscopy in one case demonstrated abnormalities in the collagen fibril organization and structure. The genomes of the two unrelated affected cats and one of the affected siblings were sequenced and individually compared to 54 feline control genomes. We searched for private protein changing variants in known human EDS candidate genes and identified three independent heterozygous COL5A1 variants. COL5A1 is a well-characterized candidate gene for classical EDS. It encodes the proα1 chain of type V collagen, which is needed for correct collagen fibril formation and the integrity of the skin. The identified variants in COL5A1 are c.112_118+15del or r.spl?, c.3514A>T or p.(Lys1172*), and c.3066del or p.(Gly1023Valfs*50) for cases no. 1, 2&3, and 4, respectively. They presumably all lead to nonsense-mediated mRNA decay, which results in haploinsufficiency of COL5A1 and causes the alterations of the connective tissue. The whole genome sequencing approach used in this study enables a refinement of the diagnosis for the affected cats as classical EDS. It further illustrates the potential of such experiments as a precision medicine approach in animals with inherited diseases.

A COL5A2 In-Frame Deletion in a Chihuahua with Ehlers-Danlos Syndrome.

Ehlers-Danlos syndrome (EDS) is a group of heterogeneous, rare diseases affecting the connective tissues. The main clinical signs of EDS are skin hyperextensibility, joint hypermobility, and skin fragility. Currently, the classification of EDS in humans distinguishes 13 clinical subtypes associated with variants in 20 different genes, reflecting the heterogeneity of this set of diseases. At present, variants in three of these genes have also been identified in dogs affected by EDS. The purpose of this study was to characterize the clinical and histopathological phenotype of an EDS-affected Chihuahua and to identify the causative genetic variant for the disease. The clinical examination suggested a diagnosis of classical EDS. Skin histopathology revealed an abnormally thin dermis, which is compatible with classical EDS. Whole-genome sequencing identified a heterozygous de novo 27 bp deletion in the COL5A2 gene, COL5A2:c.3388_3414del. The in-frame deletion is predicted to remove 9 amino acids in the triple-helical region of COL5A2. The molecular analysis and identification of a likely pathogenic variant in COL5A2 confirmed the subtype as a form of classical EDS. This is the first report of a COL5A2-related EDS in a dog.

Independent DSG4 frameshift variants in cats with hair shaft dystrophy.

Investigations of hereditary phenotypes in spontaneous mutants may help to better understand the physiological functions of the altered genes. We investigated two unrelated domestic shorthair cats with bulbous swellings of the hair shafts. The clinical, histopathological, and ultrastructural features were similar to those in mice with lanceolate hair phenotype caused by loss-of-function variants in Dsg4 encoding desmoglein 4. We sequenced the genomes from both affected cats and compared the data of each affected cat to 61 control genomes. A search for private homozygous variants in the DSG4 candidate gene revealed independent frameshift variants in each case, c.76del or p.Ile26fsLeu*4 in case no. 1 and c.1777del or p.His593Thrfs*23 in case no. 2. DSG4 is a transmembrane glycoprotein located primarily in the extracellular part of desmosomes, a complex of adhesion molecules responsible for connecting the keratin intermediate filaments of neighbouring epithelial cells. Desmosomes are essential for normal hair shaft formation. Both identified DSG4 variants in the affected cats lead to premature stop codons and truncate major parts of the open-reading frame. We assume that this leads to a complete loss of DSG4 function, resulting in an incorrect formation of the desmosomes and causing the development of defective hair shafts. Together with the knowledge on the effects of DSG4 variants in other species, our data suggest that the identified DSG4 variants cause the hair shaft dystrophy. To the best of our knowledge, this study represents the first report of pathogenic DSG4 variants in domestic animals.

ABHD5 frameshift deletion in Golden Retrievers with ichthyosis.

Ichthyoses are hereditary skin disorders characterized by the formation of scales and defects in the outermost layer of the epidermis. In dogs, at least six different breed-specific ichthyoses including a relatively common PNPLA1-related autosomal recessive ichthyosis in Golden Retrievers are known. In this study, we investigated 14 Golden Retrievers with scales that were not homozygous for the mutant PNPLA1 allele suggesting a genetically distinct new form of ichthyosis. Histopathological examinations showed lamellar, orthokeratotic hyperkeratosis, and mildly hyperplastic epidermis that led to the diagnosis of a nonepidermolytic ichthyosis. Combined linkage and homozygosity mapping in 14 cases and 30 nonaffected family members delimited a critical interval of ∼12.7 Mb on chromosome 23. Whole-genome sequencing of an affected dog revealed a single protein-changing variant within this region that was not present in 795 control genomes. The identified variant is a 14 bp deletion in the ABHD5 gene (c.1006_1019del), leading to a frameshift and altering the last 14 codons p.(Asp336Serfs*6). The genotypes at this variant showed perfect cosegregation with the ichthyosis phenotype in a large family comprising 14 cases and 72 controls. ABHD5 encodes an acyltransferase required for lipid metabolism. In humans, variants in ABHD5 cause Chanarin-Dorfman syndrome, a neutral lipid storage disease with ichthyosis. Our data in dogs together with the knowledge on the effects of ABHD5 variants in humans strongly suggest ABHD5:c.1006_1019del as candidate causative genetic variant for a new canine form of ichthyosis, which we propose to designate as Golden Retriever ichthyosis type 2 (ICH2).

A Missense Variant in SLC39A4 in a Litter of Turkish Van Cats with Acrodermatitis Enteropathica.

In a litter of Turkish Van cats, three out of six kittens developed severe signs of skin disease, diarrhea, and systemic signs of stunted growth at 6 weeks of age. Massive secondary infections of the skin lesions evolved. Histopathological examinations showed a mild to moderate hyperplastic epidermis, covered by a thick layer of laminar to compact, mostly parakeratotic keratin. The dermis was infiltrated with moderate amounts of lymphocytes and plasma cells. Due to the severity of the clinical signs, one affected kitten died and the other two had to be euthanized. We sequenced the genome of one affected kitten and compared the data to 54 control genomes. A search for private variants in the two candidate genes for the observed phenotype, MKLN1 and SLC39A4, revealed a single protein-changing variant, SLC39A4:c.1057G>C or p.Gly353Arg. The solute carrier family 39 member 4 gene (SLC39A4) encodes an intestinal zinc transporter required for the uptake of dietary zinc. The variant is predicted to change a highly conserved glycine residue within the first transmembrane domain, which most likely leads to a loss of function. The genotypes of the index family showed the expected co-segregation with the phenotype and the mutant allele was absent from 173 unrelated control cats. Together with the knowledge on the effects of SLC39A4 variants in other species, these data suggest SLC39A4:c.1057G>C as candidate causative genetic variant for the phenotype in the investigated kittens. In line with the human phenotype, we propose to designate this disease acrodermatitis enteropathica (AE).

Variants Affecting the C-Terminal Tail of UNC93B1 Are Not a Common Risk Factor for Systemic Lupus Erythematosus.

Systemic lupus erythematosus (SLE) is a heterogeneous multifactorial disease. Upregulated TLR7 signaling is a known risk factor for SLE. Recently, it was shown that specific genetic variants in UNC93B1 affect the physiological regulation of TLR7 signaling and cause characteristic autoimmune phenotypes with monogenic autosomal recessive inheritance in mutant mice and dogs. We therefore hypothesized that homologous variants in the human UNC93B1 gene might be responsible for a fraction of human SLE patients. We analyzed 536 patients of the Swiss SLE Cohort Study for the presence of genetic variants affecting the C-terminal tail of UNC93B1. None of the investigated patients carried bi-allelic UNC93B1 variants that were likely to explain their SLE phenotypes. We conclude that genetic variants affecting the C-terminal tail of UNC93B1 are not a common risk factor for SLE. It cannot be excluded that such variants might contribute to other heritable autoimmune diseases.

A COL7A1 Variant in a Litter of Neonatal Basset Hounds with Dystrophic Epidermolysis Bullosa.

We investigated three neonatal Basset Hound littermates with lesions consistent with epidermolysis bullosa (EB), a group of genetic blistering diseases. A clinically normal bitch was bred to her grandfather by artificial insemination. Out of a litter of seven puppies, two affected puppies died and one was euthanized, with these puppies being submitted for diagnostic necropsy. All had multiple bullae and ulcers involving the nasal planum and paw pads, as well as sloughing claws; one puppy also had oral and esophageal ulcers. The complete genome of one affected puppy was sequenced, and 37 known EB candidate genes were assessed. We found a candidate causative variant in COL7A1, which encodes the collagen VII alpha 1 chain. The variant is a complex rearrangement involving duplication of a 107 bp region harboring a frameshift deletion of 7 bp. The variant is predicted to truncate more than 75% of the open reading frame, p.(Val677Serfs*11). Targeted genotyping of this duplication confirmed that all three affected puppies were homozygous for the duplication, whereas 12 unaffected Basset Hounds did not carry the duplication. This variant was also not seen in the genomes of more than 600 dogs of other breeds. COL7A1 variants have been identified in humans and dogs with dystrophic epidermolysis bullosa (DEB). The identified COL7A1 variant therefore most likely represents the causative variant and allows the refinement of the preliminary EB diagnosis to DEB.

Welfare Assessment in Calves Fattened According to the "Outdoor Veal Calf" Concept and in Conventional Veal Fattening Operations in Switzerland.

The "outdoor veal calf" system was developed to encounter the demand for a veal fattening system that allows for reducing antimicrobial use without impairing animal welfare. Management improvements including direct purchase, short transportation, vaccination, three-week quarantine in individual hutches, and open-air housing in small groups in a roofed, straw-bedded paddock with a group hutch were implemented in a prospective intervention study (1905 calves, 19 intervention and 19 control farms, over one year): antimicrobial use was five times lower in "outdoor veal" farms compared to control farms (p < 0.001), but it was crucial to ensure that antimicrobial treatment reduction was not associated with decreased animal welfare, i.e., that sick animals were not left untreated. Welfare was assessed monthly on the farms, and organs of 339 calves were examined after slaughter. Cough and nasal discharge were observed significantly (p ≤ 0.05) less often in intervention than in control farms, mortality (3.1% vs. 6.3%, p = 0.020) and lung lesion prevalence (26% vs. 46%, p < 0.001) were lower; no group difference was seen in abomasal lesion prevalence (65% vs. 72%). Thus, besides reduced antimicrobial use, calf health and welfare were improved in "outdoor veal calf" farms in comparison to traditional operations.

LAMB3 Missense Variant in Australian Shepherd Dogs with Junctional Epidermolysis Bullosa.

In a highly inbred Australian Shepherd litter, three of the five puppies developed widespread ulcers of the skin, footpads, and oral mucosa within the first weeks of life. Histopathological examinations demonstrated clefting of the epidermis from the underlying dermis within or just below the basement membrane, which led to a tentative diagnosis of junctional epidermolysis bullosa (JEB) with autosomal recessive inheritance. Endoscopy in one affected dog also demonstrated separation between the epithelium and underlying tissue in the gastrointestinal tract. As a result of the severity of the clinical signs, all three dogs had to be euthanized. We sequenced the genome of one affected puppy and compared the data to 73 control genomes. A search for private variants in 37 known candidate genes for skin fragility phenotypes revealed a single protein-changing variant, LAMB3:c.1174T>C, or p.Cys392Arg. The variant was predicted to change a conserved cysteine in the laminin ß3 subunit of the heterotrimeric laminin-322, which mediates the binding of the epidermal basement membrane to the underlying dermis. Loss-of-function variants in the human LAMB3 gene lead to recessive forms of JEB. We confirmed the expected co-segregation of the genotypes in the Australian Shepherd family. The mutant allele was homozygous in two genotyped cases and heterozygous in three non-affected close relatives. It was not found in 242 other controls from the Australian Shepherd breed, nor in more than 600 other controls. These data suggest that LAMB3:c.1174T>C represents the causative variant. To the best of our knowledge, this study represents the first report of a LAMB3-related JEB in domestic animals.

Novel Brown Coat Color (Cocoa) in French Bulldogs Results from a Nonsense Variant in HPS3.

Brown or chocolate coat color in many mammalian species is frequently due to variants at the B locus or TYRP1 gene. In dogs, five different TYRP1 loss-of-function alleles have been described, which explain the vast majority of dogs with brown coat color. Recently, breeders and genetic testing laboratories identified brown French Bulldogs that did not carry any of the known mutant TYRP1 alleles. We sequenced the genome of a TYRP1+/+ brown French Bulldog and compared the data to 655 other canine genomes. A search for private variants revealed a nonsense variant in HPS3, c.2420G>A or p.(Trp807*). The brown dog was homozygous for the mutant allele at this variant. The HPS3 gene encodes a protein required for the correct biogenesis of lysosome-related organelles, including melanosomes. Variants in the human HPS3 gene cause Hermansky-Pudlak syndrome 3, which involves a mild form of oculocutaneous albinism and prolonged bleeding time. A variant in the murine Hps3 gene causes brown coat color in the cocoa mouse mutant. We genotyped a cohort of 373 French Bulldogs and found a strong association of the homozygous mutant HPS3 genotype with the brown coat color. The genotype-phenotype association and the comprehensive knowledge on HPS3 function from other species strongly suggests that HPS3:c.2420G>A is the causative variant for the observed brown coat color in French Bulldogs. In order to clearly distinguish HPS3-related from the TYRP1-related brown coat color, and in line with the murine nomenclature, we propose to designate this dog phenotype as "cocoa", and the mutant allele as HPS3co.